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Abstract 

Recent observation of large amplitude Langmuir waveforms during a Type III event in the solar 
wind have been interpreted as the signature of the electrostatic decay of beam-driven Langmuir waves. 
This mechanism is thought to be a first step to explain the generation of Type III radio emission. The 
threshold for this parametric instability in typical solar wind condition has been investigated through 
ID- IV Vlasov-Poisson simulations. We show that the amplitude of the observed Langmuir beat-like 
waveforms is of the order of the effective threshold computed from the simulations. The expected 
level of associated ion acoustic density fluctuations have also been computed for comparison with 
observations. 



1 Introduction 

During a solar flare, high energy electrons (1- 
100 keV) are expelled from the solar corona and 
travel along the interplanetary magnetic field lines, 
producing a bump on the local electron distri- 
bution function. Langmuir waves then grow via 
the so-called bump-on-tail instability. Langmuir 
waves are observed at amplitude large enough 
{E'^/{8TrnT) ~ lO'^ - IQ-") to further excite elec- 
tromagnetic waves through non-linear processes. 
These electromagnetic waves are the main signa- 
ture of Type III radio bursts. Wave- wave inter- 
action through parametric instabilities have been 
shown to be the physical mechanism underlying the 
non linear evolution of large amplitude Langmuir 
waves. 

Langmuir electrostatic decay (hereafter LED) - 
also called Langmuir decay instability (LDI) or 
parametric decay instability (PDI) in the litera- 
ture - enables energy transfer from a finite ampli- 
tude Langmuir wave L toward a second Langmuir 
wave i' and an ion acoustic wave (hereafter lAW) 
S through a three-wave resonant interaction: 

L^L' + S 



the generation of solar wind type III radio emis- 
sion at twice the plasma frequency (T2/ as a 
result of a coalescence of the two Langmuir waves 
Ginzhurg and Zheleznyakov {195i)\ : 



L 



L' 



Spectral observations of high frequency in- 
tense Langmuir waves and low-frequency ion 
acoustic waves during Type III bursts have 
been interpreted as a signature of Lang- 
muir electrostatic decay \Lin et aL\ . Further- 
more, waveform observations during Type III 
bursts of modulated Langmuir wave pack- 
ets on one hand ICairns and Robinson{1992) 



Hospodarsky and Gurnett{l995) Bale et al. (1996) 
Li et g/. (2003)1 and of lAW frequency drift asso- 
ciated with Type III electron beam velocity drift 
on the other hand Cairns and Robinson{1995) 
have both been interpreted by the Langmui 
electrostatic decay. 

The LED is a resonant parametric instability. To 
be resonant, the three- wave coupling requires the 
fundamental equations of energy and momentum 
conservation to be satisfied: 



k,, +ks 



(1) 



This process is thought to be a first step toward where uj and k are the frequency and wavenumber 
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of the waves. Moreover, for the LED to develop, 
the electric field of the mother Langmuir wave 
has to reach a critical value. The analytical 
calculation of this threshold has been discussed in 
Nishikawa{1968)\ 



Sagdeev and Galeev{1969) 



Dysthe and Franklin{l970) 



Bardwell and Gold'man(\97&)\ 
Robinson et a/.(1993)| by considering three 
monochromatic waves. The underlying idea is that 
for the instability to develop, the growth rate j^^d 
of the two product waves must be higher than 
their own linear Landau damping rates 7 , and 73 : 



The threshold for parametric decay of the 
mother Langmuir wave is thus evaluated in 
term of the electric to kinetic energy ratio 
Bardwell and Goldman{1976)\ : 



with ujl', the angular frequency of the daugh- 
ter waves. The threshold for electrostatic decay has 
been estimated for typical solar wind parameters to 
be {eoE^)/{nkBT) > 2.5 xlQ-^ [Lin et a/.(19"86b)] . 
To our knowledge, only LED that couples three 
monochromatic waves have been considered ana- 
lytically. Indeed, the analytical treatment of reso- 
nance becomes complicated when considering a fi- 
nite frequency bandwidth, in which case a numeri- 
cal treatment is necessary. 

Earlier related works on simulations of LED 
have been based on the Zakharov equations 
Sprague and Jej/er(1995j| 



Gibson et aL(1995) 



Soucek et a/. (2003)} ^ considering the instability 
as a fiuid mechanism. Recently, kinetic simula- 
tions using PIC codes [Matsukiyo et aL(2004) 
Huang and Huang{2008)\ and Vlasov codes 



Goldman et oL(1996) Umeda and /to(2008)| have 
shown that the beam-plasma interaction can satu- 
rate through LED and have been able to reproduce 
the modulated shape of Langmuir waves. 

If the threshold for LED is commonly thought to 
be at moderate amplitude, a recent numerical work 

Umeda and /te(2008)| claims that no LED would 



occur until the electric energy is comparable to the 
plasma kinetic energy. As a result, it would be no 
more possible to consider the LED as a mechanism 
for the generation of Type III radio emission. 



Henri et aL(2009)] , hereafter Paper I, recently 
reported direct observations of Langmuir waves de- 
caying into secondary Langmuir waves and acous- 
tic waves during a Type III solar event, from 
STEREO/WAVES data. They found that the 
Doppler-shifted frequencies of the three observed 
waves satisfy the resonant relations of momen- 
tum and energy conservation expected for three- 
wave coupling. A bicoherence analysis confirmed 
the phase coherence of the three waves. In this 
former work, the LED threshold and the growth 
rate of lAW density fluctuations generated by LED 
were both evaluated from analytical solutions in- 
volving a purely monochromatic three-wave cou- 
pling Sagdeev and Galeev {1969)\ . However, obser- 
vations show that: (i) the large amplitude Lang- 
muir waves are isolated wave packets with a packet 
width of the order of a few wavelengths; (ii) ion and 
electron temperatures are close (0.1 < Te/Ti < 10) 
so that ion acoustic waves associated with the LED 
should be Landau-damped. Thus several questions 
remain open about the generation and the dynam- 
ics of the LED observed in the solar wind: 

1. What are the threshold and the growth rate of 
the LED when electron and ion temperatures 
are close? What is the effective threshold of 
LED when the mother Langmuir wave is a lo- 
calized wave packet? 

2. What is the saturation level for I AW den- 
sity fluctuations in these conditions? What is 
the expected level of lAW density fluctuations 
where saturation is not observed? 

The goal of this paper is to answer these 
questions by studying the dynamics of the LED 
through ID-IV Vlasov-Poisson simulations. 

The paper is organized as follow. The Vlasov- 
Poisson simulation model is described in section [21 
The numerical results are presented in section |21 
first for a monochromatic mother Langmuir wave 
f sect. [XT]) ■ then for a mother Langmuir wave packet 
(sect. 13. 2p . Growth rates, saturation levels for 
lAW density fluctuations and threshold for LED 
are studied. The simulation results are finally dis- 
cussed in section |4] in the context of solar wind 
observations. Particular attention is paid to the 
case of equal temperature for electrons and ions, in 
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which case the growth of lAW should be hmited by 
its Landau damping. 



2 Full Vlasov-Poisson simula- 
tion model 

In typical solar wind conditions, the ratio between 
electron and ion temperatures fluctuates around 
1. The I AW produced during three- wave coupling 
should then be suppressed by Landau damping. 
This would limit the development of the lAW 
and thus the LED. However, the lAW Landau 
damping rate in such temperature condition 
cannot be solved analytically, so that it cannot 
be included in a fluid code. Kinetic effects due 
to wave-particle interactions are to be taken into 
account self-consistenly as a possible limitation for 
the instability to grow. A Vlasov-Poisson approach 
has thus been used to study LED dynamics. It 
allows to consider self-consistently the decay of the 
Langmuir wave together with damping effect on 
the product waves. Since solar wind electrons are 
weakly magnetized {uicel^pe — 10~^), magnetic 
effects are discarded. 

The Vlasov-Poisson system is solved for the elec- 
tron and ion distribution function, fe{x,v,t) and 
fi{x,u,t), with the numerical scheme described in 
[Mangeney et aL(2002)| , limiting our study to the 
ID-IV case. The equations are normalized by 
using the following characteristic electron quan- 
tities: the charge e, the electron mass nie, the 
electron density rig, the plasma (angular) fre- 
quency LUpe — Y/47m^^£?/m^, the Debye length 
Xd ^ ■\/Te/47rriee^, the electron thermal velocity 
vth.e — ^D^pe — \/Te/m^ and an electric flcld 
E = meVth,e'^pe/ Then, the dimensionless equa- 
tions for each species read: 



dt 



,dfe 
dx 



dv 



= 



(2) 



+u^ + -[E + E^:r)^ = (3) 



dt dx fj, 



du 



f^d.^ f.du;E^-f^ (4) 



where v (resp. u) is the electron (resp. ion) ve- 
locity normalized to the electron thermal velocity. 
Here /i = mi/rrie — 1836 is the ion-to-electron mass 
ratio, (p and E are the self-consistent electric po- 
tential and electric field generated by the plasma 
charge density fluctuations according to Poisson 
equation (Eq. 15). and E^"}"'' are 'external' 

drivers added to the Vlasov equation that can be 
switched on or off during the run. The electron 
(resp. ion) distribution function is discretized in 
space for < a; < Lx, with = 5000 Ad the 
total box length, with a resolution of dx — Xd- 
The electron velocity grid ranges over —5 vth,e ^ 
V < +5 vth,e, with a resolution of dv — 0.04 Vfh,e- 
(resp. —5 uth,i < u < +5 Uth,i, with a resolution 
of du = 0.04 Uth,i for the ion velocity grid, where 
Uth,i is the ion thermal velocity). Finally, periodic 
boundary conditions are used in the spatial direc- 
tion. 

In all the runs the following initial conditions 
have been taken: electron and ion distributions 
functions are initially Maxwellian with respect to 
velocity, with a random noise in density: 



fe{x,v) 



/27r 



1 + e'^^cos{kx + ^k)\ (5) 



e'^cos(fcx + V^)] (6) 



where Q — Ti/T,, is the ion-to-electron tempera- 
ture ratio set to 6 = 0.1 or 1 in the different runs. 
V'fe and 'i/'fe are random phases with a uniform 
distribution. The parameter e (resp. e') is the 
amplitude of the initial electron (resp. ion) density 
level of noise. The parameters (e,e') are chosen so 
that the development of the instability happens 
relatively rapidly, as Vlasov codes have a very 
low level of numerical noise. The influence of the 
level of initial density fluctuation on the LED 
has been studied for values 10^^ < e, e' < 10"*^. 
Neither the growth rate of density fluctuations, 
nor its saturation level are modified by the choice 
of the parameters e and e', as long as they remain 
weak. Only the time needed for the instabil- 
ity to saturate is modified. The chosen values 
e = e' = 10^^ are a good compromise to limit the 
computation time. However, when the electron 
and proton temperatures are of the same order 
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the initial perturbation on ion density is rapidly 
damped out and the generation of the lAW starts 
from the numerical noise. When starting from a 
monochromatic Langmuir wave (section 13.11) , the 
interaction time between the waves is 'infinite'. In 
this case the lAW can grow from the numerical 
noise to significant values providing we wait for 
long enough. On the other hand, when starting 
from a Langmuir wave packet (section I3.2p the 
interaction time between the waves is now finite. 
A continuous injection of noise is then needed 
to seed the instability. This is why an external 
driver E]!;°l^'^ is added in the Vlasov equation for 
ions (Eq. . The aim is to control the generation 
of a continuous incoherent noise in the proton 
density. For self-consistency, the driver is used in 
both cases = 0.1 and = 1. The details of 
the forcing are described in appendix [Xj From a 
physical point of view, lAW grow from the proton 
density fluctuations background. Note that such 
density irregularities are observed in the solar wind 
Celnikier et al.{1987)\\Harvey et aL(1988)| . 



In space conditions, the full physical process 
for the nonlinear evolution of large amplitude 
Langmuir waves is characterized by two successive 
steps. The first one is the generation of Lang- 
muir waves from the bump-on-tail instability, 
the second one the electrostatic decay of these 
Langmuir waves if their amplitudes reach the LED 
threshold. Other authors have already studied 
the generation of beam driven Langmuir wave 



packets \Omura et aZ.(1994)| Omura et aL(1996) 


Silin et a/. (2007) 




J/meda (2007) . Among their 


results, they 


lave shown that the beam- 



driven Langmuir waves are localized pack- 
ets. This localization is explained through 
the mechanism of kinetic localization in the 
framework of a nonlinear trapping theory 



Muschietti et a/.(1995) 
Akimoto et a/.(1996)j 



Muschietti et aL(1996) 
Usui et al.{2Qm)\ . It isTas 



also been shown that the beam-driven Langmuir 
wave packets are formed on time scales consider- 
ably shorter than those of parametric instabilities 
Intrator et al.{l9M)\\Akimoto et qL(1996)| . 



Our aim here is to focus the attention on the LED 
process. Motivated by the fact that the generation 
and localization of beam-driven Langmuir waves 
are decoupled from the LED process, we choose 
not to generate Langmuir waves by a bump-on-tail 



instability, but to resonantly grow the initial Lang- 
muir wave by means of an electric field E™f added 
in Eq. [51 This method enables to have a direct con- 
trol on the energy and the spectrum of the initial 
Langmuir wave and avoid other effects due to the 
non-Maxwellian character of the initial distribution 
function. Switched on at the beginning of the run, 
it acts as a driver to resonantly grow the wave with 
the desired spectrum and electric field amplitude. 
Details on the external driver are given in 

Appendix |B] The external driver is then switched 
off and the generated Langmuir wave evolves self- 
consistenly. In section [5TT1 a monochromatic Lang- 
muir wave of wavelength and amplitude is 
generated by imposing E""^ — eI,1\ (cf. Eq. [Tf)) . 
In section 13.21 a Langmuir wave packet of mean 
wavelength , packet width A and maximum am- 



init 
ext 



E. 



(2) 



plitude Ej^ is generated by imposing E[ 
(cf. Eq. [IHD. The external electric field is 
switched off when the amplitude (resp. maximum 
amplitude) of the generated Langmuir wave (resp. 
wave packet) reaches the desired value Ej^. This 
happens typically for time t < 300 ui~^ , small com- 
pared to the decay time scale. So, the resonant gen- 
eration of the Langmuir wave (resp. wave packet) 
does not interfere with the LED mechanism. 

In order to compare the simulation results with 
observations of electric waveforms during Type III 
events, the Langmuir wavelength and amplitude 
are set as indicated by solar wind observations. 
Langmuir waves grow by resonance with an elec- 
tron beam at phase velocity — uj^/k^^ ~ Vbeam ■ 
Typical electron beams span in the range 0.05 — 
0.2 c \Dulk et aL(l"987)| \Hoang et aL(1994)] . Tak- 
ing into account the Langmuir dispersion relation, 
we deduce Langmuir wavelengths in the range A^ = 
[100 — 600] \d- However, for a temperature ratio 
= 1 (resp. = 0.1), the daughter Langmuir 
wave packet is expected, from Eq. [T] and respec- 
tive dispersion relations, to be backscattered for 
\^ < 285 Xd (resp. A^^ < 385 A_d), and scattered 
forward for A^ > 285 A_d (resp. A^ > 385 A_d). 
Thereby, mother and daughter Langmuir Doppler- 
shifted frequencies are very close and hardly sepa- 
rated on observation. Thus, to be able to compare 
the simulations with observations, we choose Lang- 
muir wavelengths that allow observations of LED, 
namely A^^ = 100 A^i for the initial monochromatic 
Langmuir wave in section 13.11 and 50 A^i < A^^ < 
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400 Ad for the initial Langmuir wave packet in 
section 13.21 Since the phase velocities of the two 
expected Langmuir waves are large enough, their 
Landau damping can be neglected. Thus the phase 
velocity of the Langmuir waves do not need to be 
resolved in the electron distribution function veloc- 
ity box. On the other hand, the expected lAW 
travel at the ion-sound speed Cs = \J {T^ + Ti)/mi, 
which is of the order of the ion thermal speed when 
electron and ion temperature are close, thus leading 
to a large Landau damping. The ion-sound speed 
is resolved in the ion distribution function velocity 
box. 

The typical Langmuir waves amplitude is directly 
given by waveform observations and normalized to 
the electron temperature: 

y nkbTe 

In Paper I, the typical observed values are 10~^ < 
< 10^^. We choose several amplitudes in the 
range 10"'' < < 1. 

3 Numerical results 

During the first part of the simulation (0 < i < to) 
the external electric pump generates a Langmuir 
monochromatic wave (resp. wave packet) at the de- 
sired amplitude. At time t = tg, the external pump 
is switch off when the Langmuir wave (resp. wave 
packet) reaches the desired value. For t > to the 
Langmuir wave evolves self-consistently. We con- 
sider time t = to as the beginning of the numerical 
experiment. 

The typical intensity of the electric field of the 
lAW generated by LED is much below the level of 
Langmuir electric noise. In the following, the ion 
density is therefore used as a tracer for lAW. 

3.1 Electrostatic decay of a 
monochromatic Langmuir wave 

In this section, we consider the evolution of a 
monochromatic wave generated by the external 
pump defined in Eq. 1171 Table [T] summarizes 
the expected wavelength, wavenumbers, phase 
velocity and group velocity in the case of three 
monochromatic resonant waves, as a function 



of the mother Langmuir wavelength A^^ and the 
temperature ratio Q. 



A broad spectrum of daughter waves generated 
by the instability is clearly seen on Fig. [T] It shows 
the electric field and ion density spectrum, black 
and red lines respectively, at given times. Dashed 
vertical lines indicate the expected wavenumbers 
of the LED products as reported in table [T] In top 
left panel, t = 10'^ Wpe^, the spectrum corresponds 
to the "initial condition" where only the Langmuir 
wave at kj^ — 0.063 A^^ and its harmonic at 
0.13 A^^ are present. Then, the growth of the 
LED produces: (i) a daughter Langmuir waves L' 
at fc^, ~ 0.04 A^^ in the electric field spectrum; 
(ii) an lAW at ~ 0.10 A^^ in the ion density 
spectrum. Both are shown in the next two panels, 
i = 5 X 10'' uj-^ (top right panel) and t = 10^ cj"/ 
(bottom left panel). Finally, at t > 10^ (bottom 
right panel) harmonics of the lAW are produced, 
\ow-k fluctuations are generated and the LED 
saturates. We need to stress here that during 
the decay phase, even if the Langmuir wave L is 
monochromatic, the product waves L' and S are 
both wave packets with wavenumbers centered 
on the expected wavenumber. This results from 
the fact that different fc-channels are available for 
energy transfer from the Langmuir wave toward 
its decay products. Indeed Fig. [5] shows the linear 
growth rate vs. k for different wavenumbers of 
the daughter Langmuir wave electric field (in 
black) and lAW density fluctuations (in red) 
spectrum. Growth rates for both daughter waves 
are overplotted with fc^-axis and fc^,-axis such 
that kg + k^, = fc^. As expected, the growth 
rates are the same for each couple of product 
waves that verify the wavenumber resonant 
condition. The spread in wavenumber is about 
Afc^, — Akg ~ 0.025 for both waves, thereby a 
relative spectral spread of Ak^/k^ ~ 0.2 for the 
lAW and Afc^, /k^, ~ 0.6 for the Langmuir product 
wave. Note that the generation of a large spectrum 
of product waves limits the spatial coherence of the 
interaction. This point is discussed in section [4.21 



We then compute an 'integrated' (over space) ion 
density fluctuations < S n > generated by LED 
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Waves 


Wavelength 


Wavenumber 


Phase velocity 


Group velocity 






(Ad) 




(wth = Adw-I) 


i^th = Ar>w-i) 


e = 0.1 


L 


= 100 


= 0.063 


= 16 


vl ^ 0.189 




L' 


A^, = 132 


k^, = (-)0.047 




= -0.141 




S 


As 57 


fcg = 0.110 


Cs = 0.024 


= 0.024 


e = 1 


L 


A^ = 100 


k^ = 0.063 


= 16 


vl = 0.189 




L' 


A^, = 154 


k^, = (-)0.041 


w^, = -24.5 


= -0.123 




S 


A, = 60 


fcg = 0.104 


0.033 


Cs = 0.033 



Table 1: Expected waves involved in the Langmuir electrostatic decay, given the initial Langmuir wave- 
length and using equations [T] for 9 = 0.1 (top line) and = 1 (bottom line). It shows for each wave 
respectively the expected wavelength (in Debye length Ad), corresponding wavenumber (in inverse Debye 
length Xj^) used in the paper figures, phase and group velocities (in A/jwZ.^) 



defined by: 



<5n> (t) = W 



{n{t,x) -no)^dx (8) 



with the mean ion density no = 1 in dimensionless 
units. We define the saturation level 6nsat as 
the maximum value reached by the mean density 
fluctuations < Sn > (t) during the simulation. The 
growth rate for average density fluctuation is here 
iLBn = 2.3 X lO-^Wpe, or 7^^„ = 4.1 x lO'^/s 
when expressed in term of the lAW frequency, 
defined by = c/Ag = 5.5 x 10"^. This 
means that the characteristic time scale for the 
density to grow is about 24 lAW periods, in this 
case. The saturation level Susat = 8.7 x 10~^, 
obtained at time t = 1.4 x 10^ ^pe, is of the 
order of the expected saturation level defined by 
the ratio of electric energy to the thermal energy 
5n°,,^^El/{T, + T,) = 2.3xm-\ 

The simulation has been repeated for several ini- 
tial Langmuir wave amplitude, with values in the 
range 10"^ < < 1, and for two temperature 
ratios 6 = 0.1 and 6 = 1. 

We found that the evolution becomes strongly 
non-linear when the Langmuir wave electric field is 
"not small", i.e. when > 0.3 for 9 = 1 and 
Ej^ > 0.2 for 8 = 0.1. When strong non linear pro- 
cesses take place the evolution is not simply driven 
by the LED mechanism. Conversely, for small val- 
ues of E^ , the system is mainly driven by the LED 
process so that we can compute the growth rate for 
average density fluctuations and the density fluctu- 
ations saturation level as previously deflned. 



The growth rate deduced from simulations 7^ed 
of the average density fluctuations < Sn > deduced 
from the simulation is displayed in Fig.[3]for the two 
temperature ratios. As shown before, for too high 
values of the Langmuir amplitude - about Ej^ ~ 
0.3, corresponding to an electric to thermal energy 
ratio of 0.1 - strong nonlinear effects arise before 
LED. Also, as expected, the growth rates are lower 
for 8 = 1 than 8 = 0.1, due to the increase of 
lAW Landau damping in the first case. Contrary 
to former results Umeda and Jto(2008)| , we show 
here that the LED is observed for initial Langmuir 
waves with electric energy five orders of magnitude 
lower than the plasma thermal energy. 

An analytical effective growth rate 7 
has been overplotted (dashed lines) 
fined as the difference between the analytical 
full monochromatic case of undamped monochro- 
matic waves in an homogeneous background 7*')jj-, 
Sagdeev and Galeev{196^ and the lAW Landau 



analytical 

LED 

It is de- 



damping 7^°"''""- 



analytical 

'LED 



' LED 



, Landau 



where 



„th 



nksT me 



■ip_\i/i 



(9) 



(10) 



7 is of the order of the lAW frequency with so- 
lar wind parameters. The growth rate computed 
from Vlasov simulations ^^^d is one to two orders 
of magnitude lower than ^<™^aiyt^cai _ r^j^^ discus- 
sion concerning the discrepancy between the ana- 
lytical monochromatic case and Vlasov simulations 
is postponed to section [ 
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r 


a 




Tp/Te = 0.1 


0.026 


1.11 


0.59 


Tp/Te = 1 


0.025 


1.82 


0.30 



Table 2: Numerical values for the fit for the growth 
rate of LED driven density fluctuations in expres- 
sion [TTJ 

We also looked for the dependence of the growth 
rate with respect to the mother Langmuir wavevec- 
tor kj^ . Finally, the growth rate for density fluctua- 
tions generated by LED has been fitted by a power 
law: 

^,,,=TE'^X (11) 
The fitting parameters F, a and j3 are given in 
table m 

The saturation level of density fluctuation Suaat 
is summarized in figure |4] for = 0.1 (top panel) 
and = 1 (top panel). The expected saturation 
level < (^ng"* >, expressed as the initial Langmuir 
electric energy to the total kinetic energy ratio, 

<<Sng-*>=ii;2/(T, + T0 (12) 

is overplotted (dashed lines). The obtained 
saturation level of density fluctuation is in good 
agreement with the expected values in both cases. 

To summarize, the simulations of the Langmuir 
Electrostatic Decay from an initial monochromatic 
Langmuir wave have shown that: 

1. the threshold for the instability to grow, ex- 
pressed in term of the Langmuir wave elec- 
tric energy, is at least 5 orders of magnitude 
lower than the plasma thermal energy when 
O.K e < 1; 

2. the product waves are resonantly generated 
over a broad range of wavenumbers, naturally 
producing narrow wave packets; 

3. growth rates of lAW density fluctuations are 
one to two orders of magnitude lower than 
the analytical values deduced from the pure 
monochromatic case; 

4. saturation levels for lAW relative density fluc- 
tuations are of the order of the ratio of Lang- 
muir electric energy to the total kinetic energy. 



3.2 Electrostatic decay of a Lang- 
muir wave packet 

In the following, we consider the evolution of 
a finite-amplitude Langmuir wave packet, gener- 
ated by the external electric field pump defined in 
Eq. [THl We recall that, in order to mimic the pres- 
ence of a low level of small scale ion acoustic tur- 
bulence, we introduce a source term of ion acous- 
tic noise. This forcing generates incoherent proton 
density fluctuations at a level 5n/n ~ 10~^ much 
smaller than the expected level of density fluctua- 
tions generated by LED. 

Figure [S] shows the LED of a Langmuir wave 
packet with wavelengths centered on = 200 A^, 
a packet width A = 2000 A^i, a maximum initial 
electric field E^^ — 6 x 10~^, and the two temper- 
ature ratio = 0.1 (left panels) and = 1 (right 
panels). The mother and daughter Langmuir 
wave packets can be followed in the top panels 
that show the space-time evolution of the electric 
field density energy £'(x,t)^/2. The Langmuir 
mother wave packet propagates towards the right 
and emits backscattered Langmuir wave packets. 
The bottom panels shows the temporal evolution 
of ion density fluctuations during the decay of 
the Langmuir wave. lAW density fluctuations are 
generated locally (by ponderomotive force from the 
two Langmuir wave packets beats) and propagate 
forward at the ion sound speed (dashed line). 
When = 0.1 lAW propagate and escape the area 
where LED develops. Conversely, when = 1, 
lAW density fluctuations are heavily damped as 
soon as the waves escape the area where LED 
occurs. Thereby, the waves can be observed only 
where the mother Langmuir wave decays. Finally, 
as for the case of monochromatic Langmuir waves, 
the LED growth rate is lower for a larger value of 0. 

The simulation has been repeated for different 
values of and mean wavelength A^^ , with a 
packet width of A = 10 A^^ . Figure [6] summarizes 
the results for the evolution of the Langmuir wave 
packet, each point representing a single simulation. 

During the decay process, the mother wave 
packet generates a daughter wave packet travelling 
at a different group velocity. The region of pa- 
rameters leading to electrostatic decay is displayed 
with green squares in flgure [SI Eventually the two 
wave packets separate thus stopping the LED pro- 
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cess. Therefore, LED is efficient only if the inter- 
action time between the two wave packets is longer 
than the growth time for the daughter waves. The 
growth time of the instability is controlled by the 
Langmuir electric field amplitude: the larger the 
amplitude, the smaller the growth time (Fig. 
The interaction time is controlled by the Langmuir 
wavelength: the larger the wavelength, the smaller 
its group velocity and so the larger the interaction 
time. Thus, for low amplitude and/or short wave- 
length, the Langmuir wave packet propagates at 
its group velocity without non-linear interactions. 
The region of parameters leading to a linear be- 
havior without electrostatic decay is displayed with 
orange crosses in figure El The efficiency of LED 
(green) region increases toward the linear (orange) 
part (i.e. LED is observed for lower A^^ and E'^ ) if 
temperature ratio Q decreases or/and the size A of 
the mother wave packet increases. Finally, strong 
non linear effects dominate the evolution of high 
amplitude and/or large wavelength Langmuir wave 
packets, and LED is no more the dominant process. 
The region of parameters where other non linear 
effects are dominant is shown with red crosses in 
figure El These strong non linear effects seem to be 
the signature of strong turbulence: Langmuir col- 
lapse and formation of cavitons. This evolution is 
out of scope of the present study and will be studied 
in a future work. 

The effective threshold for the electrostatic de- 
cay of a Langmuir wave packet can be estimated 
by imposing that the interaction time Tint is equal 
to the inverse growth rate of daughter wavepack- 
ets 7^£,o . The growth rate for LED 7^£,^ has been 
obtained from simulations of monochromatic Lang- 
muir waves LED and fitted in Eq. [TTl The interac- 
tion time Tint is evaluated from: 

nnt^A/{vl~vl,) (13) 

where (w^ — v^, ) is the difference of group velocity 
between the two waves that first separate ~ here the 
mother and daughter Langmuir waves - and A the 
packet width of the mother Langmuir wave, which 
is assumed to be about the length of interaction. 
Using Eq. [T] and ~ 2 kj^, Eq. [T31 can be written 
in normalized units as: 

T„t~A/(6fcJ. (14) 

Finally, Tint = I^^d gi'^^s the effective LED thresh- 
old of a Langmuir wavepacket with wavevector fc^ 



and a packet width A: 

rpthreshold _ I 
LED " V AF / ^ ' 

expressed in normalized units. The semi analyti- 
cal threshold is overplotted in Fig. [51 (red line) in 
order to validate its dependency with the mother 
Langmuir wave vector kj^ . 

We also performed a series of simulations in or- 
der to validate the dependency of E^^^^'^^°^'^ with 
A. This time the simulation starts with a mother 
Langmuir wave packet of mean wavelength Al — 
200 but with different packet width values in the 
range 5Al < A < 30 Al, covering the typical 
range for Langmuir wavepackets in the solar wind, 
and initial amplitude E^. The results are presented 
in Fig. [71 together with the semi analytical thresh- 
old (red line). In both cases, Eq. [TSl is in agree- 
ment with simulations of LED of localized Lang- 
muir wave packets. 

The threshold decreases when: (i) the ion to elec- 
tron temperature ratio decreases, since the Lan- 
dau damping of lAW decreases and their effec- 
tive growth rate increase; (ii) the mean Langmuir 
wavelength increases, since Langmuir wave pack- 
ets with longer wavelengths propagate at smaller 
group velocities thus increasing the available inter- 
action time with the ion background; (iii) the width 
of the Langmuir wave packets increases, since the 
interaction time between the mother and daughter 
waves increases. The effective threshold obtained 
by Vlasov simulations and described in this section 
is compared to observations of LED in the solar 
wind in the next section. 

Finally, in two recent works of ID- 
IV Vlasov-simulations, [£/meda(2007)| and 
\Umeda and /to(2008)| reported the evolution 
of beam-excited Langmuir waves. Their simu- 
lation parameters are Q = 0.1, a beam speed 
Vbeam = 8 X Vth,e with a beam temperature equals 
to the electron core temperature and a beam 
density ratio of 0.1% and 0.5% respectively. These 
beams generate Langmuir waves at a phase velocity 
of about Vbeam - Vth,e = 7 X Vth,e, 1-6. a Laugmuir 
wavelength Al — 40 x Xd- The amplitude of the 
Langmuir waves reaches E^ ~ 0.3 in the first case, 
El ~ 1 in the second case. They observed the 
LED process in the second case only. We have 
overplotted their results in Fig. El (black cross and 
diamond respectively). These previous simulations 
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agrees with Eq. [151 for these parameters, the 
equation predicts that LED of locahzed Langmuir 
wave packets should occur for high Langmuir 
electric energy levels (of the order of unity) . 

To summarize the simulations of the Langmuir 
electrostatic decay from an initial Langmuir wave 
packet: 

1. we have illustrated how the localization of 
Langmuir wave packets is crucial for the evo- 
lution of Langmuir decay, by limiting the in- 
teraction time between mother and daughter 
waves; 

2. we have shown that the ion acoustic waves are 
generated locally, where the wave packets in- 
teract; then lAW are damped as soon as they 
escape the region where the mother and daugh- 
ter Langmuir waves interact when electron and 
ion temperature are equal, or escape this re- 
gion when the electron temperature is higher 
than the ion temperature; 

3. we used results from the Langmuir electro- 
static decay of a monochromatic Langmuir 
wave fsection IO]) to compute a semi analytical 
threshold for the electrostatic decay of Lang- 
muir wave packets (Eq. llSp ; this semi analyt- 
ical threshold has been shown to be in agree- 
ment with simulations of the electrostatic de- 
cay of Langmuir wave packets. 



4 Discussions 

In the following, we discuss, first, the initialization 
of LED in the case = 1, second, the discrepancy 
between the LED growth rate obtained from sim- 
ulations and analytical estimation. LED threshold 
obtained from simulations is then compared to ob- 
servations. Finally, we show that the saturation of 
the instability gives an upper limit to the expected 
level of observed density fluctuations. 

4.1 Initialization of Langmuir elec- 
trostatic decay 

In the case of equal electron and ion temperature, 
associated with a strong Landau damping of ion 



acoustic fluctuations, one may wonder whether pro- 
cesses complementary to the resonant interaction of 
the waves could facilitate the LED. In particular, a 
decrease (even local) of the lAW Landau damping 
during the beginning of the LED process could ease 
the initialization of the instability. Such decrease 
can happen either through two main processes, (i) 
Wave-particle interactions due to the trapping of 
ions in the lAW potential well could modify the 
ion distribution function by forming a plateau at 
the ion sound speed which, in turn, could then de- 
crease the Landau damping rate. However, no such 
plateauing is observed, (ii) The beats of the mother 
and daughter Langmuir waves, could heat the elec- 
trons so that the ion-to-electron temperature ratio 
decrease locally, leading to a partial suppression 
of the lAW Landau damping. However the tem- 
perature ratio during the simulation has too small 
variations A9/0 < 1% to really modify the lAW 
Landau damping. This hypothesis is thus also ruled 
out. 

Since no such complementary processes are 
present in the simulations, the initiation of the LED 
is likely to be caused by the resonant interaction 
of the waves, that dominates locally the Landau 
damping of lAW. 

4.2 Growth rate for Langmuir elec- 
trostatic decay 

We have seen in sect. 13.11 that a large discrepancy 
exists between the growth rate of mean density fluc- 
tuations obtained from the simulations, 7^^^^ j ^^.d 
the analytical one deduced from three monochro- 
Sagdeev and Gofee^ (1969)0 



matic waves. 



th 

' LED 



The existence of a strong lAW Landau damping 
7i,„„.„„, such that 7^^^ = 7*'^^ - 7i,„„.„„^ could 
explain this difference. Since Landau damping for 
lAW cannot be found analytically when electron 
and temperature are close, we performed comple- 
mentary simulations and measured an lAW Landau 
damping of ^^^^^ 
(resp. 7i,„„da„ = 
much lower than ^th 

ILED 'led iLandc 



-1.9 X 10"^ ujpe for e = 0.1 
8.9 X 10^-* LUpe for 8 = 1), 
Thus, the hypothesis that 
is to be ruled out. 



^The influence of the resolution in velocity has been 
checked, and the results are unchanged as far as the resolu- 
tion in velocity is not too low. The results are also indepen- 
dent of the size of the box, as long as the box remains larger 
than the coherence length of the daughter wave packets. 
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Note that the growth rate measured from sim- 
ulations ah-eady takes into account the effect of 
the Landau damping. Thus, the difference between 
^LED ^"^d ryojudyticai giiould be Carried by the eval- 
uation of 7* Actually, the main difference be- 
tween the analytical treatment of three-wave reso- 
nance and Vlasov simulations is that the daughter 
waves are treated as monochromatic in the first case 
whereas they are observed to have a non negligible 
spectral width in the second (see Fig [1] and [2]). Sim- 
ulations have shown that the spectral width of the 
product waves increases when: (i) ~ 1, in which 
case the lAW dispersion relation is numerically ob- 
served to spread out around the analytical branch; 
therefore the resonance can occur with (w, fc)- values 
slightly different from the theoretical expectation; 
(ii) the energy of the Langmuir mother wave in- 
creases, in which case the linear approximation for 
a (5-shaped resonance is no longer valid. The limited 
spatial coherence of daughter wave packets implies 
that their growth is localized where mother and 
daughter wave packets interact, thus strongly limit- 
ing their growth rate. This explains why the growth 
rate deduced from simulations is much lower than 
the analytical one. 

4.3 Threshold of Langmuir electro- 
static decay and type III obser- 
vations. 

The TDS observation mode of the S/WAVES 
experiment on board the STEREO mission 

Bougeret et aL(2007)| gives access to in-situ elec- 
tric field waveform in 3D with an equivalent spec- 
tral resolution up to 60 kHz. In Paper I, we have 
shown evidence for nonlinear coupling between 
Langmuir waves at about 10 kHz and ion acoustic 
waves at about 0.2 kHz. We recall here the global 
plasma parameters for these observations. The 
one-hour- average wind speed from WIND/3-DP 

Lin et aL(1995)| is about Vsw = 315 km s^^. 
The electron temperature observed by WIND/3DP 
is Te ~ 10^ K , and the electron density in the 
solar wind, estimated from the plasma frequency, 
is about Tie — 10^ m~^. From the electron density 
and temperature, the Debye length is A^i ~ 20 m. 

Figure |8] displays the observed value of Langmuir 
electric field, normalized as described in Eq. [71 for 
the whole data set of waveforms where LED has 



been observed. Their wavelength, not directly mea- 
sured, is evaluated as follow. By taking into ac- 
count the Doppler-shift caused by the solar wind 
and noting that the lAW should propagate here 
anti-sunward, the ion acoustic wavelength reads: 

\ \ ^ \ I f observed 

As = {Vsw + Cs)/ Js 

where the sound speed Cs is evaluated by 

Cs = \JkB{Te +Tp)/mp, 

where ks is the Boltzman constant, rup the proton 
mass and Tp the proton temperature that remains 
unknown but is of the order of Tg in the solar wind. 
The mother Langmuir wavelength A^^ is then about 
twice the LED-produced ion acoustic wavelength. 

However, caution should be taken when directly 
comparing this threshold value with the observa- 
tion. Indeed solar wind type HI Langmuir wave 
packets are most probably localized in 3D. The ob- 
served waveforms are ID spatial cuts of the real 
3D structures and thus only give a lower limit on 
the width of the wave packet, as well as its max- 
imum amplitude. We assume that the Langmuir 
wave packets have a 3D gaussian shape, and that 
the spacecraft crosses it somewhere within a dis- 
tance from the center of the 3D structure of the 
order of the full width at half maximum FWHM. 
Therefore a realistic range for the maximum electric 
field El'^'"'- reached by the 3D wave packet is calcu- 
lated as: Ef" < El"""^ < e Ef", where iJf " is the 
observed level of Langmuir electric field. Data are 
plotted with such error bars. 

The LED threshold expressed in Eq. [15] is over- 
plotted for two values of the ion to electron tem- 
perature ratio Q — 0.1 and Q = I, using computed 
values from table [51 The observed electric field am- 
plitudes are of the order of the threshold computed 
from Vlasov-Poisson simulations - in the case of a 
localized wave packet and « 1 - confirming that 
the development of LED is compatible with the ob- 
served events. 

4.4 Saturation of LED: expected 
level of observed density fluctu- 
ations 

In our simulations, for initial Langmuir wave pack- 
ets in the range of amplitudes that correspond to 
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observations of Langmuir waves during type III 
events, the lAW-like density fluctuations associated 
with the LED do not reach the saturation level. In- 
deed, the mother and Langmuir wave packets that 
propagate at different group velocities separate be- 
fore the saturation stage, thus stopping the growth 
of the lAW. This is confirmed in the observations, 
since most of the observed waveforms show: (i) 
mother and daughter Langmuir waves with differ- 
ent energies; (ii) lAW with no harmonics. The sat- 
uration level < 5ng"* > /no = El/{Te + Ti) should 
thus be an upper limit on the expected level of 
density fluctuations on observations of LED dur- 
ing type III. 

Density fluctuations at the frequencies consid- 
ered here can be measured in space from the vari- 
ations of the floating potential of the spacecraft 
\Peders e w(1995)] when it crosses the region where 
LED is observed. The voltage signal observed on 
STEREO /Waves waveforms with monopole anten- 
nas channels contains information on the floating 
potential of the spacecraft. This remains however 
to be calibrated. Knowing the saturation level of 
lAW during LED, and given the calibration, we 
could check whether the observed level of density 
fluctuations is consistent with our simulations. 

5 Conclusion 

In order to study the origin of electromagnetic ra- 
dio emissions during type III bursts, we have re- 
ported in this paper ID-IV Vlasov-Poisson simu- 
lations of the Langmuir Electrostatic Decay. The 
simulations have been done in typical solar wind 
conditions: ratio of the electron to ion temperature 
from 0.1 to 1, mother Langmuir wavelengths typ- 
ical of those observed during type III events and, 
most important, by considering localized Langmuir 
wave packets. The main results are the following: 

1. Langmuir electrostatic decay develops even 
when the electron and ion temperatures are 
close. Its threshold, when considering a 
monochromatic wave, is at least flve orders of 
magnitude lower than the plasma thermal en- 
ergy when 0.1 < < 1. 

2. Langmuir electrostatic decay resonantly gen- 
erates daughter waves over a broad range 
of wavenumbers, naturally leading to narrow 



wave packets. This limits the length of coher- 
ence which is why the growth rate is one to 
two orders of magnitude lower than the ana- 
lytical values deduced from a pure monochro- 
matic case. 

3. The behavior of daughter ion acoustic waves 
depends on the temperature ratio. Ion acous- 
tic waves can escape the region where the 
resonant coupling takes place and propagate 
when electron temperature is higher than pro- 
ton temperature. Conversely, they are damped 
as soon as they escape the resonant coupling 
area when temperatures are equal. 

4. We conflrm that the saturation level for lAW 
density fluctuations is of the order of the ratio 
of Langmuir electric energy to the total kinetic 
energy. However, for the range of amplitudes 
that correspond to observations, the lAW-like 
density fluctuations associated with the LED 
should not reach the saturation level because 
(i) the mother and Langmuir wave packets that 
propagate at different group velocities separate 
before the saturation stage, (ii) harmonics of 
the lAW are seen at saturation in the simula- 
tions but, to our knowledge, not in Type III 
observations. 



Finally, an effective threshold has been ob- 
tained (Eq. lisp for localized Langmuir wave 
packets and compared to STEREO/WAVES 
observations. The observed Langmuir 
electric field during type III reported in 
Henri et aZ.(2009)| is in the range of LED 



effective threshold computed from Vlasov sim- 
ulations, thus confirming the interpretation of 
these observed electric field waveforms in term 
of the LED of type III beam-driven Langmuir 
waves. 

The physical mechanism responsible for the 
generation of electromagnetic radio waves during 
type III burst is still under study. The process de- 
scribed by Ginzburg and Zheleznyakov {1958)] and 
leading to the generation of type III radio emis- 
sion at twice the local plasma frequency is a two 
step process, (i) First the beam-driven Langmuir 
wave decays through LED, (ii) then the mother and 
daughter Langmuir waves coalesce to generated the 
electromagnetic wave. 
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Up to now, observations have shown that res- 
onant couphng between Langmuir waves and ion 
acoustic waves does occur during a type III burst 
Henri et a/. (2009)] . The present paper also shows 



that Langmuir electrostatic decay may occur in 
solar wind conditions and that the threshold is 
reached by the observed Langmuir electric ampli- 
tude. These two complementary studies thus con- 
firms step (i) does occur. 

Langmuir Electrostatic Decay generates two 
counter-propagating Langmuir waves in opposite 
direction, but step (ii) requires obliquely propagat- 
ing Langmuir waves for waves to couple and pro- 
duce the transverse electromagnetic wave. Does 
the density inhomogeneities scatter the Langmuir 
waves enough to introduce a significant perpendicu- 
lar component to their wave vector? An important 
step would be to study the interaction of Langmuir 
waves with an inhomogenous background, in order 
to check this hypothesis. Finally, the coalescence 
of counter-propagating Langmuir waves is hard to 
observe because the signal would be hidden by the 
generation of the harmonic of the mother Langmuir 
wave. Simulations of step (ii) with input from the 
observations could at least show whether the coa- 
lescence process could occur in the solar wind. 



Appendix: Generation 
the ion density noise 



of 



The continuous injection of ion density noise is 
driven by an external fields E"°l'^'^ added to the 
ion dynamics (Eq|3]). This driver is defined by: 



5t = 500 xti;"^'^ (about 80 plasma oscillations). This 
means that for an lAW of wavelength Xia = 100, 
the forcing lasts 1/20*'' of a period, much shorter 
than the typical lAW oscillation time in order to 
generate an incoherent noise. This forcing thus gen- 
erates density fiuctuations much smaller than the 
level of density fluctuations generated by LED in 
our simulations. 



B Appendix: Generation of 
the initial Langmuir wave 



init 
ext ' 



In this appendix, we describe the driver 
added in the Vlasov equation for electrons (Eq. [5]) . 
This 'external' electric field controls the generation 
of the initial Langmuir wave. It acts as a driver 
that resonantly grow a Langmuir wave propagat- 
ing in only one direction, with the desired spec- 
trum and electric field amplitude. In section 13.11 
a monochromatic Langmuir wave of wavelength A^^ 
and amplitude Ej^ is resonantly excited by the ex- 



ternal electric field pump by defining E'^ 



init 
ext 



E 



(1). 



Eil\ (x, t) = E^^^ (t) cos{kox - Loot) (17) 



In section 13.21 a Langmuir wave packet of mean 
wavelength , packet width A and maximum am- 
plitude E^ is resonantly excited by defining El^f — 



E. 



(2). 



(2) 

E^J^{x^i) — Eo{t) cos{kQX-'UJot) exp 



In both cases, the pump wavevector k[ 







(18) 
27r/A, 



^ext i^J — ^ext 



lons^vaax 



cos{2'KX / X)cos(uj\t 



I cos(2ttx/ X)cos{uj\t + ^ 
(16) 

It introduces a fiat spectrum for wavelength over 
the range 50 < A < 1000. E'^°l^<'{x,t) is nor- 
malized in order to have a maximum amplitude 
of j^^J'J^^'"^"-^ = 1 X 10^^. The frequencies uj\ are 
chosen to satisfy the dispersion relation of lAW, 
u!\ = (27r/A)cs with Cg = a/(1 -I- 0)//x the ion 
sound speed. The phases ip'il{t) have a step-like 
variation, constant over a time interval St. At the 
end of each interval, they are independently drawn 
from a uniform distribution. This way, the genera- 
tion of ion acoustic noise is made of a succession of 
time coherent forcing for time intervals of duration 



^p a^^Y^ the frequency wq = y^l + 3fc§ are chosen to 



TOsfy the Langmuir dispersion relation, in order 
TO generate a Langmuir wave at the desired wave- 
length Xj^. The width of the wave packet is set 
to A = 10 A^. The pump wave packet moves at 
the Langmuir group velocity — Sfco/wo from its 
initial position xq = L^/S. The time-dependant 
external electric field amplitude Eq (t) scales on the 
desired Langmuir ampitudc E^^ : 



E^,'\t)^r^E,{l-exp{-t/T)). 



(19) 



In the monochromatic case (resp. wave packet 
case), the parameter r is set to 30 lu~^ (resp. 
100 LOpJ-) so that the external electric field ampli- 
tude smoothly increases during a few Langmuir 
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oscillations, in order to avoid the generation of 
other plasma waves induced by step-like switch of 
the pump. The parameter ry is set experimentally 
to 5 X 10^^ (rcsp. 10^^), so that the characteristic 
time scale to resonantly grow the Langmuir wave 
(resp. wave packet) is large compared to its 
oscillation time scale, but small compared to the 
decay time scale. So, the resonant generation of 
the Langmuir wave (resp. wave packet) does not 
interfere with the LED mechanism. 
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Figure 1: Spectrum of both electric field (in black) and ion density (in red) at 4 different times (indicated 
on each panel), from simulation of LED of a monochromatic Langmuir "wave "with -wavenumber — 0.063 
and initial electric amplitude Ej^ — 3 x 10~^. Dashed lines indicate the expected wavenumbers of the 
LED "waves. Top left panel displays the initial conditions and the bottom right panel sho-ws the spectrum 
at saturation. Note the presence of the second harmonic of the lAW in this last case. 




Figure 2: Spectral growth rate for the Langmuir "wave L' electric field (in red, fc^,-axis at the bottom, 
left-ward) and the lAW density fluctuations (in black, /c^-axis at the top, right-ward) generated by the 
LED of a monochromatic Langmuir -wave {E^ = 3 x 10^^ and 9 = 1). Both axis are scaled in order 
to show the resonant condition in wavenumbers: fc^, + k^, ~ k^ . Vertical error bars are 1-cr error bars 
for the fitting of the growth rate. The vertical dotted line indicates the expected wavenumbers for 
monochromatic LED products. 
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Figure 3: Growth rate for lAW density fluctuations vs. Langmuir wave initial electric field amplitude 
for two temperature ratios. Top panel: 6 = 0.1; Bottom panel: 9 = 1. The growth rate is ex- 
pressed in plasma frequency unit uipe (left axis) and in JAW frequency fj^ (right axis). The full line 
is a fit of Vlasov simulations (stars). The dashed line shows the analytical growth rate ^^^a'y'^cai from 
Sagdeev and Galeev{196^ when LED-generated lAW remain monochromatic including the linear Lan- 
dau damping of lAW. The vertical dotted line sets the limits above which strong nonlinear effects occurs 
before LED starts. 
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Figure 4: Average density fluctuation at saturation Susat (expressed in % of mean density) for different 

Langmuir wave initial electric field amplitude, for both = 0.1 (top panel) and = 1 (bottom panel). 
The dashed lines show the expected level of saturation feg"* = (0.5 x E-^)/{Te + Ti) defined as the 
Langmuir electric energy to the total kinetic energy ratio. 
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Figure 5: Top panels: space-time evolution of the electric energy density; Bottom panels: space-time 
evolution of the ion density fluctuations. Left panels: Q = 0.1, right panels: = 1. Full lines show the 
expected group velocities of the Langmuir wave packets while dashed lines displays the ion sound speed. 




Figure 6: Evolution of a Langmuir wave packet, with different wavelengths centered on A^^ (bottom 
axis) but a same packet width of A = 10 A^^, and different initial electric fields i?^ (left axis). Left 
Panel: 9 = 0.1; Right panel: 6 = 1. In green: simulations where LED is observed; In orange: 
simulations where the available interaction time between the waves is lower than the LED timescale; In 
red: simulations where strong NL effect arise before/instead of LED. Previous simulation results form 
\Umeda{2007) Umeda and /to(2008)| are overplotted in the left panel with black cross and diamond 
respectively. The red line is the semi analytical threshold expressed in Eg [T5] 
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Figure 7: Evolution of a Langmuir wave packet, with a mean wavelength = 200 but different packet 
widths A (bottom axis, normalized to lambda j^), and different initial electric fields E^^ (left axis). Left 
Panel: O = 0.1; Right panel: = 1. In green: simulations where LED is observed; In orange: 
simulations where the available interaction time between the waves is lower than the LED timescale; In 
red: simulations where strong NL effect arise before/instead of LED. The red line is the semi analytical 
threshold expressed in Eg [T5] 
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Figure 8: Observed Langmuir electric field amplitude - normalized like the simulations - vs. wavelength. 
The threshold for LED computed from Vlasov simulations is overplottcd for two values of the temperature 
ratio e = 0.1 and 9 = 1. 
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